Undulatory swimming in viscoelastic fluids 
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The effects of fluid elasticity on the swimming behavior of the nematode Caenorhabditis elegans 
are experimentally investigated by tracking the nematode's motion and measuring the correspond- 
ing velocity fields. We find that fluid elasticity hinders self-propulsion. Compared to Newtonian 
solutions, fluid elasticity leads to up to 35% slower propulsion speed. Furthermore, self-propulsion 
decreases as elastic stresses grow in magnitude in the fluid. This decrease in self-propulsion in 
viscoelastic fluids is related to the stretching of flexible molecules near hyperbolic points in the flow. 

PACS numbers: 47.63.Gd, 47.15.G-, 47.20.Gv, 83.50.Jf 



Many microorganisms have evolved within complex 
fluids, including soil, intestinal fluid, and human mu- 
cus [lHaj- The material properties or rheology of such 
fluids can strongly affect an organism's swimming be- 
havior. For example, in the case of freely swimming 
spermatozoa, the flagcllum shows a regular sinusoidal 
beating pattern. Once the organism encounters a vis- 
coelastic medium, this regular beating pattern is replaced 
by high-amplitude, asymmetric bending of the flagellum. 
The motility behavior of the sperm cell is affected by 
its fluidic environment 0, 0|, which in turn can affect 
human fertility A major challenge is to understand 
the mechanism of propulsion in media that displays both 
solid- and fluid-like behavior, such as viscoelastic fluids. 

Our current understanding of swimming at low 
Reynolds (Re) numbers is derived mainly from investi- 
gations in Newtonian fluids [8-12[. Here Re = p UL/p, 
where p and p are the fluid density and viscosity, and U 
and L are the organism's speed and characteristic length 
scale. At low Re, locomotion results from non-reciprocal 
deformations in order to break time-reversal symmetry; 
this is the "scallop theorem" It has been recently 

shown 0, E| that the scallop theorem may break down 
for viscoelastic fluids due to the fluid's history-dependent 
stresses that grow nonlinearly with strain rate. These 
elastic stresses can dramatically change the flow behav- 
ior even at low Re [14j . 

The effects of fluid elasticity on swimming at low Re 
have been considered in theory and numerical simulation. 
For an infinite waving sheet immersed in a second-order 
fluid flBj . it was shown that elasticity augments propul- 
sion speed. Recently, it was shown that for the case of 
an infinite undulating sheet Q and cylinder Q , viscoelas- 
ticity decreases swimming speed compared to Stokesian 
Newtonian cases. By contrast, a two-dimensional nu- 
merical simulation for a finite undulating sheet using the 
Oldroyd-B model [l6j showed that fluid elasticity could 
in fact augment swimming speed when the beating fre- 
quency / is equal to the inverse of the fluid relaxation 
time A; that is, the Deborah number De = fX = 1. 
Despite these recent efforts, there is a dearth of experi- 



mental investigations of swimming in viscoelastic fluids, 
and the effects of fluid elasticity on swimming are still 
not clear. 
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FIG. 1. Color Online, (a) Sample snapshot of the nema- 
tode C. elegans swimming in buffer solution (p — 1.0 mPa-s). 
Lines represent nematode's "skeleton" and its centroid path, 
(b) Corresponding contour plots of the nematode's bending 
curvature k(s, t) over 3 swimming cycles, (c) Streamlines 
computed from instantaneous velocity fields of Newtonian 
(Re < 10~ 3 ) and (d) polymeric (Re < 10~ 3 ;De = 3.0) flu- 
ids. Arrows in (c,d) indicate flow direction and the box in (d) 
shows a hyperbolic point in the flow. 

In this letter, the effects of fluid elasticity on an 
undulatory swimmer are experimentally investigated at 
low Re by tracking the swimmer and tracer particles 
in the flow (Fig. 1). The organism is the nematode 
Caenorhabditis elemns, a roundworm widely used for bi- 
ological research [171 ] that swims by generating traveling 
waves [III, 18 1 . Overall, we find that fluid elasticity hin- 
ders propulsion compared Newtonian fluids (Fig. 2) due 
to the enhanced resistance to flow near hyperbolic points 
for viscoelastic fluids. 

Experiments arc performed in fluid-filled channels that 
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are 15 mm wide and 600 /im deep. The swimming motion 
of C. elegans is imaged using bright-field microscopy and 
a CMOS camera at 125 frames per second. The nema- 
tode is approximately 1 mm long and 80 /im in diameter. 
The objective focal plane is set on the longitudinal axis 
of the nematode body. All data presented here pertain 
to nematodes swimming at the center plane of the fluidic 
channel. Out-of-plane recordings are discarded. An av- 
erage of 15 nematodes is recorded for each experiment. 
More detailed information on the experimental methods 
can be found in [lj| [20[. Figure 1(a) shows a sample 
snapshot of a nematode swimming in a water-like solu- 
tion at Re = 0.2 as well as the nematode's shape-line or 
"skeleton" and its centroid path. Here, swimming speed 
U is calculated by differentiating the nematode centroid 
position over time. 

Newtonian fluids of different shear viscosities [i are pre- 
pared by mixing two low molecular weight oils (Halo- 
carbon oil, Sigma- Aldrich). Viscoelastic fluids are pre- 
pared by adding small amounts of carboxymethyl cellu- 
lose (CMC, 7 x 10 5 MW) into water 0. CMC is a 
long, flexible polymer with an overlap concentration of 
approximately 10 4 ppm . The polymer concentration in 
solution ranges from 1000 to 6000 ppm resulting in fluid 
relaxation times A that range from 0.4 s to 5.6 s, re- 
spectively. These solutions are dilute and do not show 
significant shear-thinning viscosity [lj|, particularly in 
the range of typical swimming shear-rates of 1 to 20 s _1 . 
Nevertheless, in order to rule out the effects of shear- 
rate dependent viscosity, an aqueous solution of the stiff 
polymer Xanthan Gum (XG) that is shear-thinning but 
possesses negligible elasticity is also used in experiments. 

An important quantity that is used to characterize the 
swimming behavior of undulatory swimmers, such as C. 
elegans, is the bending curvature defined as n(s,t) = 
d<fi/ds. Here, is the angle made by the tangent to 
the x-axis in the laboratory frame at each point along 
the body centerline, and s is the arc length coordinate 
spanning the head of the nematode (s = 0) to its tail 
(s = L). Figure 1(b) shows the spatio-temporal evolu- 
tion of the nematode's body curvature k(s, t) for 3T, or 3 
swimming cycles. The contour plots show the existence 
of periodic, well-defined diagonally oriented lines charac- 
teristic of bending waves, which propagate in time along 
the nematode body length. From such contour plots, we 
can extract kinematic quantities such as the nematode's 
swimming frequency / and wavelength X w as well as the 
wave speed c = X w f. For the nematode shown in Fig. 
lb, / ~ 2 Hz, X w = 2.5 mm, and c = 5 mm/s. 

The flow fields produced by the swimming nematode 
are investigated using particle tracking velocimetry poj ]. 
Examples of streamlines computed from instantaneous 
velocity fields are shown in Fig. 1(c) and 1(d) for the 
Newtonian and viscoelastic cases, respectively. Here, 
Re < 10 -3 for both fluids, and De = 3.0 for the vis- 
coelastic fluid. Overall, the streamlines display large re- 
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FIG. 2. Color Online, (a) Swimming speed, (b) bending wave 
speed, and (c) kinematic efficiency of C. elegans in Newtonian 
(red circle) and viscoelastic (blue square) fluids as a function 
of fluid viscosity. Triangle symbol represents the non-elastic 
Xanthan Gum solution. The data shows that fluid elastic- 
ity decreases the nematode's swimming speed and efficiency 
when compared to a Newtonian fluid of same viscosity. For 
fj, > 30 mPa-s, the nematode's swimming speed decreases in- 
dicating a limit in power for this type of organism. 



circulation flow structures, or vortices, that are attached 
to the nematode's body. Such patterns are similar to 
the flow visualizations of Gray and Lissmann plj who 
associated such recirculation zones with regions of max- 
imum transverse (nematode) body displacement. The 
streamlines for the Newtonian and viscoelastic cases are 
qualitatively different, with the appearance of a distinct 
hyperbolic point near the nematode for the latter case. 

We now address the question of whether fluid elastic- 
ity hampers or enhances swimming speed. The nema- 
tode's swimming speed as a function of fluid viscosity for 
both Newtonian and polymeric fluids is shown in Fig. 
2(a). For relatively low viscosity values, the swimming 
speed is independent of fluid viscosity, and the values of 
U are nearly identical for both cases. For /x > 30 mPa-s, 
the swimming speed decreases with increasing fi even for 
Newtonian fluids. The decrease in U with increasing fi at 
low Re is most likely due to the nematode's finite power. 
We note that, for a nematode swimming with constant 
power at low Re, P ~ F^ rag lJ ~ [iU 2 where P is power 
and Fdrag is the drag force the fluid is exerting on the 
nematode. Results show that, over the limited range of 
U, the C. elegans' 1 propulsion speed shows a decay that 
is slower than pT 1 ! 2 , which strongly suggests that the 
nematode does not swim with constant power. The max- 
imum power generated by the organism is approximately 
200 pW (p = 30 mPa-s) 0. 
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Nevertheless, we find that the values of U for viscoelas- 
tic fluids can be 35% lower than the Newtonian fluid of 
same shear viscosity. For example, the nematode's swim- 
ming speeds for the viscoelastic and Newtonian cases are 
0.18 mm/s and 0.25 mm/s, respectively, even though the 
viscosity for both fluids is 300 mPa-s (Fig. 2a). The de- 
crease in swimming speed in CMC (polymeric) solutions 
is not due to shear-thinning effects since nematode swim- 
ming in the non-elastic, shear-thinning fluid (XG) showed 
no apparent decrease in propulsion speed (Fig. 2a, trian- 
gle symbol) compared to the Newtonian case. 

The nematode's swimming behavior is further investi- 
gated by measuring the ben ding curvature k along the 
nematode's body centerline 12, [lij]. In Fig. 2(b), we 
show the bending wave speed c as a function of viscos- 
ity. Results show that viscoclasticity has negligible effect 
on the nematode's swimming kinematics. That is, the 
changes in kinematics including the decrease in beating 
frequency and wave speed are due to viscous effects only. 
In addition, there is no evidence of change in motility 
gait (e.g. swimming to crawling) as fi increases since the 
beating amplitudes remain constant (A 0.26 mm) even 
for the most viscous fluid (fj, = 400 mPa-s). 

Figure 2(c) shows the nematode's swimming efficiency 
as a function of fluid viscosity for both the Newtonian and 
polymeric fluids. Here, swimming efficiency is defined as 
the ratio of the swimming speed U to the bending wave 
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speed c [21[. For the Newtonian case, the swimming effi- 
ciency increases with /i until a finite asymptotic value is 
eventually approached. For CMC (polymeric) fluids, the 
efficiency initially follows the trend of Newtonian fluids 
because the fluid elastic stresses are very small (De « 0). 
At fi rj 30 mPa-s, we observe a new branch in which 
efficiency decreases with fluid viscosity. This viscoelas- 
tic branch is observed at De « 1, where the undulation 
frequency of the swimmer might couple to the fluid relax- 
ation time. Overall, the kinematic swimming data show 
that fluid elasticity hinders both the organism's swim- 
ming speed and swimming efficiency at low Re. 

The effects of fluid elasticity on the nematode's swim- 
ming behavior are best illustrated by plotting the nor- 
malized swimming speed U/Un as a function of the Deb- 
orah number (De = /A), where C//v is the Newtonian 
speed. Figure 3 shows that the normalized swimming 
speed decreases monotonically with De, and reaches an 
asymptotic value of w 0.4 as De is further increased. 
In other words, as the clastic stresses increase in mag- 
nitude in the fluid, it introduces a larger resistance to 
propulsion, therefore decreasing the nematode's swim- 
ming speed. A similar trend is observed in gels using a 
"two- fluid" model [22 1. 

Next, the experimental results on swimming speed are 
compared to recent theoretical predictions [j, |23j . We 
note that for all the experiments presented here, the ra- 
tio of the solvent viscosity to the total solution viscosity 
S = Hsoiv/^soi is below 0.05, where [i S oiv = 1-0 mPa-s is 
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FIG. 3. Color Online. Swimming speed normalized by Newto- 
nian speed (Un) as a function of Deborah number. The data 
(squares) show that propulsion speed decreases as elasticity in 
the fluid increases. Solid line shows general trend from pl.l23l| 
where 8 = 0.05 (see text). Dashed line corresponds to predic- 
tions of ||, H^] using kinematic data from this work. Dotted 
line corresponds to numerical simulations of [l(|. 



the solvent (buffer) viscosity and /i so / is the solution vis- 
cosity. For the case of an infinitely long, two dimensional 
waving sheet Q and cylinder [23| with prescribed beat- 
ing pattern, it is predicted that the swimming speed de- 
creases with increasing De. While the experimental data 
supports the predicted trend, there is still quantitative 
discrepancies between the experimental and theoretical 
results as shown in Fig. 3. Some of the possible reasons 
for the observed discrepancies may be the finite length 
of the swimmer and the assumption of small beating am- 
plitude in the theoretical works. That is, only small de- 
flections arc considered for both the waving sheet and 
cylinder while the nematode shows significant bending. 

We also compared the experimental results to a re- 
cent two-dimensional numerical simulation of a finite, 
large-amplitude waving sheet using the Stokcs-Oldroyd- 
B model [lij . The simulation predicts an interesting en- 
hancement of the sheet swimming speed at De = 1 (Fig. 
3). The experimental results do not reveal such swim- 
ming speed enhancement (Fig. 3) in viscoelastic fluids. 
For De > 1, the simulation predicts a gradual decrease 
in U. The discrepancies between the experiment and the 
simulations are most likely due to the difference in the 
swimming beating patterns. While simulations used a 
left-moving traveling wave with an amplitude that in- 
creased from head to tail, our experiments with C. ele- 
gans reveal a traveling wave with an exponential decay 
from head to tail [l8|] . 

In order to gain further insight into the effects of fluid 
elasticity on swimming, we investigate the flow fields gen- 
erated by C. elegans at Re < 10~ 3 for both Newtonian 
and viscoelastic fluids (Fig. lc,d); De = 3.0 for the vis- 
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FIG. 4. Color Online, (a) Velocity decay normal to the C. 
elegans swimming direction (inset) for Newtonian and vis- 
coelastic cases (De = 3.0) at Re < 1CP 3 . The fluid viscosity 
is n — 200 mPa-s. Lines corresponds to exponential fits (see 
text), (b) Mean square displacements of the Newtonian and 
viscoelastic fluids after 20 beating cycles T. 



coelastic case. In particular, we are interested in the 
velocity decay normal to the nematode's swimming di- 
rection because of its relevance to, for example, hydro- 
dynamic interactions and collective swimming. The ef- 
fective viscosity of both fluids is /i = 200 mPa-s (c.f. 
Fig. 2a). Figure 4(a) shows the normalized velocity mag- 
nitude of the fluid flow |V|/|V| maa; as a function of the 
normalized distance r/L away from the nematode in the 
normal direction, as shown in the inset. Here, r is the 
distance normal to the nematode with origin at the fluid- 
body interface and L is the nematode body length (1 
mm) . The velocity decays quite rapidly in less than a 
half-body length [20j, and it follows a seemingly expo- 
nential decay of the form |V"|/|V| ma a: = exp(— ^ £) pre- 
viously obtained by Lighthill for an undulating sheet Q . 
By comparison, the viscoelastic case shows a velocity de- 
cay rate that is initially faster (a = 0.56 ± 0.03) than 
the Newtonian case (a = 0.74 ± 0.04), indicating that 
elasticity hinders fluid transport around the nematode. 

The flow transport properties are further investigated 
by computing the mean square displacement (MSD) 
(Ar 2 ) of fluid particles advected in the flow for up to 
20 swimming cycles (Fig. 4b). The slope of the MSD as 
a function of time is a relative measure of particle trans- 
port due to flow. Since the Peclet number is large, 
O(10 6 ), both Newtonian and viscoelastic fluids have a 
slope k that is well above unity. That is, the fluid trans- 
port induced by the nematode swimming is non-diffusive. 
Elasticity, however, hinders fluid transport as shown by 
the lower value of k in Fig. 4(b). This is mostly likely due 
to the sudden increase of elastic stresses near regions of 
high velocity gradients such as hyperbolic points. Near 
such regions, the extensional viscosity of a solution of 
flexible polymers can be orders of magnitude larger than 
a Newtonian fluid (24[. Polymer molecules can be eas- 
ily aligned and stretched, which results in an increase 
in hydrodynamic drag along the molecules and poses an 
additional resistance to fluid transport and swimming. 



In conclusion, we have experimentally investigated the 
effects of fluid elasticity of the swimming dynamics of 
undulatory swimmers at low Re. We find that fluid elas- 
ticity hinders the propulsion of the nematode C. elegans 
at low Re. The swimming speed decreases as fluid elas- 
ticity is increased. This trend is qualitatively similar to 
theoretical and numerical results [H, 0, [l6| . Furthermore, 
elastic stresses in the fluid can alter the flow field gener- 
ated by nematodes, and the presence of hyperbolic points 
in viscoelastic flows can result in large extensional viscosi- 
ties and resistance to flow. This implies that foraging, 
feeding, and mixing may become difficult in strongly vis- 
coelastic media. We note that generally, elastic response 
is not limited to extensional viscosity effects; they could 
also take the form of the hoop stresses that are associated 
with circulating flows. We therefore expect the dynamics 
of swimming in viscoelastic media to depend very much 
on the type and strength of the swimming stroke. 
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